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It is the common lore to assume that knowing the equation for the probabihty distribution func- 
tion (PDF) of a stochastic model as a function of time tells the whole picture defining all other 
characteristics of the model. We show that this is not the case by comparing two exactly solvable 
models of anomalous diffusion due to geometric constraints: The comb model and the random walk 
on a random walk (RWRW). We show that though the two models have exactly the same PDFs, 
they differ in other respects, like their first passage time (FPT) distributions, their autocorrelation 
functions and their aging properties. 



One often thinks that perfect knowledge of the PDF as 
a function of time (achieved by better experimental tech- 
niques or by data mining) completely determines the un- 
derlying stochastic model. In what follows we illustrate 
that this is not the case by revisiting two exactly solvable 
models often invoked when explaining anomalous diffu- 
sion in labyrinthine structures of percolation type: the 
comb model introduced in 1 (which mimics trapping 
in the dangling ends) and the RWRW introduced in [2 
(which mimics the tortuosity of the chemical path). Both 
are visualized in Fig. [T] The RWRW model was intro- 
duced as a simplified (quenched) version of the single-file 
diffusion. In polymer physics, the RWRW corresponds to 
one of the regimes of reptation [3] , namely, to the motion 
of the chain as a whole along its primitive path. These 
facts make evident its relation to other models of diffu- 
sion under time-dependent constraints, which are often 
described using a fractional Brownian motion approach. 
Although both models are known for around quarter of 
a century and were quite carefully investigated, the two 
were never confronted. Doing so is however quite instruc- 
tive. Thus in the long time limit both models, the comb 
(of the CTRW class ,1 ) and the RWRW (in a class of 
itself, but a close relative of the single file diffusion) show 
exactly the same PDF which can be described by the 
corresponding fractional diffusion equation 4J. However, 
the models differ in any other respect: the comb model 
exhibits aging while the RWRW does not, the distribu- 
tion of FPTs to a given site are different and also the 
auto-correlation functions differ. 

Note that both the comb and RWRW models stem 
from Markovian random walk (RW) models on the corre- 
sponding geometries, i.e., the probabilities of steps from a 
site (i, j) on the underlying track depend only on the po- 
sition (i, j). In both cases the RWs have independent but 
not identically distributed increments. We are interested 
in the projection of the RW onto the a;-axis. The fact that 
the same value of x may correspond to different values of 
y in the comb model, or to different positions along the 
contour length n in the RWRW model (see Fig.[l]), and 
that the further motion depends explicitly on this y or n, 
introduces the memory on the previous positions. Intro- 
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FIG. 1: (color online) Left: RWRW model. The trajectory 
of a random walk (acting as a track on which a RWer may 
move) is displayed in black. Circles mark the position x = —2 
at different 'chemical distances' along the track. Right: A RW 
on a symmetric comb structure. The RW trajectory is marked 
in a lighter shade (red). The projection of the trajectory 
on the a:-axis, or backbone, leads to waiting times caused 
by the RW's sojourn in a tooth (depicted by circles of sizes 
representing the appropriate waiting time), thus a CTRW. 

ducing the memory is the cost we pay for neglecting the 
irrelevant coordinate (y or n), the dependence on which 
however does not disappear. The kind of non-Markovian 
behavior of the reduced models is strongly different: The 
diffusion on the comb falls into a semi-Markovian CTRW 
class, the one on the RW track is an anti-persistent RW, 
a close relative (but not a member of the class) of frac- 
tional Brownian motions. 

In the CTRW as represented by the comb model the 
probability density of being at position x after n time 
steps is given by the subordination expression 

poo 

= y]p(a;,n)x(n,i) ~ / p{x,n)x{n,t)dn (1) 

n=0 

where p{x, n) gives the probability of reaching x after 
n steps and xii^it) is the probability to make n steps 
within time t. Here t is measured in the units of time 
steps r = 1. Note that for t ^ r the typical number 
of steps n is large and can be considered as a continu- 
ous variable, turning the sum in Eq. ^ into an integral 
over n. The waiting time distribution along the comb 
backbone is given by the return time probability of a 
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FIG. 2: (color online) PDFs of RWRW (smooth red line) vs. 
RW on a comb (dashed blue line) for a different number of 
time steps: (a) t=100, (b) t=1000, (c) t=5000, (d) t=6500. It 
can be clearly seen that the PDFs of the two models coincide. 



RW on a tooth to the backbone, which is approximately 
V'(i) ~ (27r)-i/2i-3/2 ^ i-jjg Laplace domain, xin,t) 
and are related via 



X(n,s) = s-i[l-^(s)]^"(s), 



(2) 



where /(s) denotes the Laplace transform of a function 
f{t). Substituting V(s) ^ I - V2s'^/'^ into Eq. ^ and 
performing the inverse Laplace transform results in: 



2 / 



(3) 



a Gaussian distribution restricted to the non-negative 
half-line. Substituting Eq. ^ and p{x, n) (which tends 
to a Gaussian for large n for a regular RW on the back- 
bone) into the continuous limit of Eq. ([T]), we get: 

= / fT, exp 



X 

'2n 



exp 



n 

'2i 



dn. 



(4) 

Let us now calculate the PDF of the particles' displace- 
ment in the RWRW model, pertinent to many realiza- 
tions of walks and starting points. Let us consider a RW 
on a RW starting at a: = and consider its position after 
i 3> 1 time steps in a continuous approximation when 
the displacement n of the RWer along the trajectory of 
the RW is well-approximated by a Gaussian: 



K{n, t) 



1 



: exp 



n 
2i 



(5) 



The projection on x of this displacement along the tra- 
jectory is again given by a Gaussian depending on the 
displacement \n\: 



p{x,n) 



1 



exp 



x~ 
'2N 



(6) 



The displacement in x as a function of t is expressed by: 



pix,t) 



p{x, n)K,{n, t)dn. 



(7) 



Substituting Eqs. ^ and ^ into Eq. ([7| results in ex- 
tly the same expression as Eq. (|4| . This is numerically 
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verified in Fig. [2j The models are essentially twins! 

Let us however stress a fundamental difference between 
the two models. Eq. ^ implies that the time t is counted 
from the beginning of the walk starting at (0,0). In the 
course of time the walker typically enters deeper into the 
dangling ends. If the observation starts at a later instant 
t > 0, the walker will typically need time to return to the 
X-axis to take a step in a relevant direction. Therefore the 
waiting time PDF for the first step after t > is different 
from ^{t) d], leading to a different xin^t): the diffusion 
on the comb shows aging typical of CTRW [71 [8] . On the 
contrary, the RWRW model does not age: the same dis- 
tribution K,{n,t) holds, whenever one starts. This prop- 
erty can be inferred from Fig. [3] showing MSDs for both 
models for observations starting at different times from 
the beginning of the RW. 

Note that the movements in the a; and y directions on 
the comb are correlated: the more steps in the x direc- 
tion are made, the less time remains for steps in the y 



direction. The total MSD (r^) 



{y ) is found to 



be diffusive, {r'^{t)) ^ n, [9 , which is numerically veri- 
fied (inset in Fig.jsf. It is interesting to note that though 
the a;-projection of this RW follows a non-ergodic sub- 
diffusive process, the 2d picture corresponds to ergodic 
normal diffusion. This emphasizes the importance of vig- 
ilance on experimentalists' side when assuming that a 
projection onto a lower dimension gives a representative 
description of the process. 

Let us now turn to the properties of the process other 
than the PDF and calculate the FPT distributions for 
the two models, starting with the comb. As in a general 
CTRW, the FPT for a comb is given by the subordination 
formula 



/(i)-^/(n)^„(t) 



(8) 



71=0 



where f{n) is the FPT probability at the n-th step in a 
regular RW and ipn (t) is the probability density to make 
the n-th step at time instant t (note that the first passage 
only takes place exactly when the step is taken, so that 
?/,„(t) ^ xin, t)). The FPT of a regular RW is asymptot- 
ically given by the Smirnov distribution: 



fin) 



\/27rn3/2 



exp 



X 

2n 



(9) 



In the continuum limit one may approximate the sum by 
an integral over n. Since in the Laplace domain '0n(s) = 
V'"(s) for s — > we have tpnit) — exp(n(l — ip{s)), re- 



sulting in /(s) - f(n)t 



dn = e ^l^ll^'l , which 
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FIG. 3: (color online) MSDs in a comb and in a RWRW 
for different time lags after the beginning of the RW (r = 
0,1000,3000). The MSDs for different time lags for the 
RWRW overlap (upper, red curve). The lower (blue) curves 
represent the MSDs of the comb model. The larger the time 
lag the longer it takes the MSD to reach its asymptotic be- 
havior. The inset shows the MSD of the the 2d comb model 
averaged over time (black dashed line) and over an ensemble 
(overlapping solid blue), showing the ergodicity of the two- 
dimensional motion. The behavior of {x^) (lower solid blue) 
is shown for comparison. 



is the characteristic function of a one-sided Levy law of 
index 1/4. Its behavior ioi t ^ oo is 



fit) 



-21x1 



r(-i/4) 



-5/4 



(10) 



We stress that the same result follows from the solution 
of the fractional diffusion equation with an absorbing 
boundary [5]. 

Let us now calculate the FPT distribution for the 
RWRW model [9 . We fix the starting point x — Q and 
the final point xi. The RW track on which the RW 
then takes place has two strands connecting xq to the 
two points corresponding to xi (existing due to a recur- 
rence of RWs in Id). The distances h and I2 along the 
strands are independent random variables whose distri- 



butions p{li) and pQ.2) follow the Smirnov law, Eq. (22 1. 
Therefore the corresponding RW on the track takes place 
on a finite interval (—^1,^2) starting at Z = 0. Changing 
variables we pass to a symmetric interval (— L/2, L/2) of 
length L = I1+I2 with the starting point at y = {li — l2)/2 
inside it. The survival probability of a particle within 
this interval with absorbing boundaries (in the diffusion 
approximation, _D = 1/2) is given by 



where a„ 



00 

E 



(-l)'"a„,exp 



Dt) (11) 



4 cos 



7r{2m+l)y 



7r^(2TO+ 1)^ 

Z2 



/7r(2m -1-1). At long times 



only the first term (m = 0) matters, so that 



lcos(g)exp (-^Dt 



(12) 



where Dt = n/2. Averaging over the position of the 
starting point y and over the length of the interval gives: 

00 

m = JJ ^l,2{t)Pih)p{l2)dhdl2. (13) 





Substituting Eqs. (25) and (22 1 into Eq. i26h and chang- 



ing variables to = I / (h + I2) and z = h/h we get 



CX3 



l + z /ttI-z 



exp -TT^DtS,^ - 



4z 



Odedz. 



(14) 



Evaluating Eq.(14) for t 
where C 



0.64 



00 we get [9j: <^{t) ~ 
Ts'a number constant. Thus the survival 
probability follows the pattern of a simple Id RW with 
the FPT density 



fit) ~ ^-3/2 



(15) 



in agreement with [10], and not the slower decay of 
CTRW of f{t) (X t~^/'^. These results are numerically 
verified and shown in Fig. |4] The anomaly of diffusion 
however shows up in the x-dependence of the character- 
istic decay time tc (say the median of the survival proba- 
bility), which scales as oc x^, exactly like in the case of 
the comb, and not as t^ cx x^ like for the simple diffusion. 




FIG. 4: (color online) FPT densities for the comb (upper 
curve, blue) and for the RWRW (lower curve, red). The 
dashed lines give the theoretical asympototic behaviours. 
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Another fundamental intrinsic characteristic is the 
auto-correlation function. Let Si be the i-th step of the 
walker, so that Xn = J2l 



■11=0 



Si. The MSD is therefore: 



(16) 



i,j=0 



where the average is taken over different realizations and 



C,;, 



{siSj) is the step-step correlation function. 



Here the stationary, ergodic RWRW model differs cru- 
cially from the non-stationary CTRW one, since for sta- 
tionary RWRW processes Cij only depends on the dif- 



and is an oscillat- 



ference of its arguments Ci 
ing, slowly decaying function {vide infra). On the other 
hand, for CTRW different steps along the backbone of 
the comb are uncorrelated, so that C{i,j) = M{i)Sij 
with M{i) being the probability that the time-step t — i 
corresponds to a step on the backbone and not on a 
tooth, i.e. to the rate of relevant steps. This one can 
be evaluated via x{n, t) [5] and for large t behaves as 
M{t) ~ (27r)-i/2 ^-1/2^ r^Yve explicit dependence of 
C{i,j) on time is a fingerprint of the non-stationarity 
of CTRW. 

To asses Cr — (s„s„+t-) for a RWRW we have to look 
separately at even and odd time lags r — \i ~ j\. The 
first step s„ of the RW corresponds to zero delay and 
is taken to be in a positive direction along the track (to 
the right). Since the track of the RW is an uncorrelated 
RW itself, the two steps are correlated only if they cor- 
respond to the motion over the same step (bond) of the 
track, so that s^Sn+r = ±1 depending whether the bond 
is crossed in the same direction (at even lags) or in the 
opposite directions (at odd lags), and are uncorrelated 
otherwise. Crossing of the same bond takes place with 
probability 1/2 and only if the walker was at its starting 
position (even lags) or at its position after performing 
the first step (odd lags) . For t = 2k the probability that 
the RW on a track returns to its starting point is given 
by the binomial distribution binom(fc, 2k— 1, 1/2) (where 
binom {k, n,p) = " P)""'')- For r = 2fc + 1 

the return to the endpoint of the first step takes place 
with the probability binom(fc, 2fc, 1/2). We therefore 
have Co = 1 and Ci = —1/2, and generally: 



C2k = (l/2)binom(fc,2fc- 1,1/2) 
C2k+i = (-l/2)binom (fc, 2k, 1/2) 



(17) 



Thus a RWRW model is a bona fide anti-persistent RW 



model. Eq. (17) has been numerically verified. 



Note that both correlation functions, the one for 
RWRW and the one for the comb, lead to the same MSDs. 
For RWRW we have: (xl) = 2J2l^g{n-k)Ck-Co-Cn. 
For large t — n the two last contributions Cq — 1 and 
C„ — > can be neglected compared to the sum and 
Cfe can be approximated via the Stirling formula giving 
Cfe = ±(27rfc)~^/^. Passing from the sum to the integral 



we get for large n {x'^) — \f2Tr^v}l'^ . 



'1/2 



For the comb we 
^2/^ k-^l'^dk = 



get (x^)=ELi(2-)-^/^fc 
^2'K~^n}-l'^ , same as above. 

Let us summarize our findings. We discussed the be- 
havior of two popular models of geometrically induced 
subdiffusion: the random walk on a comb (the model of 
the CTRW class) and the random walk on the random 
walk's trajectory, belonging to the class of anti-persistent 
random walks. Although both models lead to exactly 
the same PDFs at all times, they are crucially different 
with any other respect like FPTs or correlation functions, 
which mirror the fact that the RWRW model is station- 
ary and the CTRW model is not. 

The difference in these behaviors leads us to a follow- 
ing conclusion. We say that a model is described by some 
equation if several relevant properties of the model can be 
obtained by solving this equation with corresponding ini- 
tial and boundary conditions. For example, the fractional 
diffusion equation describes the comb model because not 
only the Green's function (PDF p{x,t), the solution for 
free boundaries), but also the FPT distribution to the 
point xi follows from the corresponding equation, now 
with the boundary condition p{xi,t) — 0. This is defi- 
nitely not the case for the RWRW model, for which we 
have to admit that the knowledge of the Green's function 
(and of the equation for it) is not enough to obtain the 
FPT, and the model is not fully described by this equa- 
tion. Any conclusions about the nature of the process 
based on the PDF (or moments) alone or on the knowl- 
edge of the equation for such a PDF may thus stand on 
a quite shaky basis. 



SUPPLEMENTARY MATERIAL 
MSD of the comb model 

Let us define the total number of time steps as n and 
the number of steps taken in the x-direction as Hx . The 
number of steps taken in the y-direction is then riy = 
n — Ux- Using CTRW arguments @j and given that the 
distribution of soujourn times in the teeth is given by 
^{t) ~ ^^t~3/2^ |-]^g one- dimensional ensemble MSD 



along the backbone takes the form 



(•^ ) ens ^ ij^x) ^ Tl ^ . 



Similarly, for the y-axis we have: 

(y )ens 



n,,) = in — tIt) = n — n 



1/2 



(18) 



(19) 



The two-dimensional MSD is then the sum of the two 
one-dimensional MSDs, resulting in a linear time depen- 
dence: 



(r^) = (x^) + (y2) ^ n. 



(20) 
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This result is numerically verified, displayed in the inset 
of Fig. 3 in the main text. 

One should note that one commonly treats the motions 
in the x and y directions as independent [111 112] which 
leads to a two-dimensional propagator being the product 
of the two one-dimensional propagators. This, in turn, 
leads to a two-dimensional MSD non-linear in time, the 
one behaving as 



)ens 



{'^ }ens 



,1/2 



(21) 



Although this is asymptotically linear, the numerics in 
Fig. 3 from the main text clearly exhibit a completely 
linear time dependence at all time scales including the 
shortest ones, thus rejecting this view in favour of the 
one presented above. 



Eq. 20 describes normal diffusion on the comb struc- 
ture, which, contrary to the diffusion in the projection 
on the x-direction, is a stationary and ergodic random 
process. We also note that a random walk on a perco- 
lation cluster at criticality is subdiffusive, meaning that 
the comb structure is inherently inadequate to model a 
percolation cluster. 



The FPT distributions 

Let us obtain the FPT distribution of the RWRW 
model, again by averaging over the initial point and the 
realizations of the walk. We fix the starting point x — 
and the final point xi. The RW trajectory on which the 
random walk than takes place has two strands connect- 
ing these two points (since in ID a RW is certain to reach 
any given point on a line) , and therefore the correspond- 
ing random walk takes place on a finite interval. The 
distances li and I2 along the corresponding strands are 
independent random variables, both distributed accord- 
ing to the Smirnov law 



X 



3/2 

1:2 



exp 



4^1, 



(22) 



Therefore the corresponding RW on the RW takes place 
at an interval {—h, h) starting at Z = 0, or, equivalently, 
on an interval of length L = li + I2 starting at a point at 
a distance h from the left boundary of the interval. It 
is convenient to change the variables, and to consider a 
symmetric interval {~L/2, L/2) with the starting point 
y = {li — h)/^ inside it (being at the distance li from its 
left boundary and at the distance I2 from the right one). 
The probability to find a particle at position I within 
this interval with absorbing boundaries (in the diffusion 
approximation, D — 1/2) is given by 



2 

m=0 



where bm = TT{2m+ l)/L. Integration over I gives us the 
survival probability 



00 



$1,2 (i) = ^ (-1)"^ cos (b^y) exp i^blDt) (24) 

m— ^' 



which has to be averaged over the position of the starting 
point y and of the length of the interval. The overall ex- 
plicit evaluation of the mean of this series is tedious, but 
finding the long time asymptotic behavior is relatively 
simple. To do this, we note that at long times only the 
first term (m = 0) plays a role, so that 



— cos — exp 



L2 



Dt 



(25) 



where Dt = n/2. The average over the position of the 
starting point y and of the length of the interval is given 
by: 



00 



(26) 



We substitute Eqs. (251 and (22) in Eq. (26) and make 
a change of variables, L — li + I2 and z = I2/I1 (the 
Jacobian of this transformation is d{li,l2)/d{L, z) — 
L/(l -t- z)^). Making a further transformation, ^ — 1/L, 
leads to 



1+ Z ^TT 1 

757^cos(-- 



2 1 + 2 



d^dz. 



(27) 

We are interested in the asymptotic behavior of this 
function for i — >■ 00. The integral over ^ can be taken in 
quadratures: 



exp 



1 



2V7rDt 
and tends to 



exp 



TT^Dt^^ 

x*{l + z)^ 
Mz^-K^Dt 



x'^{l + z)' 



4z 
erfc 



a;2(l + z)2 

SZTTX^ 



2VTrDi 



(28) 



for t — > 00 (unless z gets extremely large, when it starts 
to decay as z^^). Such large values of z do not play any 
role, since the integral 



I2 



TT i 

COS ■ 



2 1 



1 + z 



dz 



(29) 

where 1F2 is 



converges: It is equal to 1F2 (^1, 5, f , - ig 
the generalized hypergeometric function (the integral can 
be transformed to the form (Eq. 2.5.8.1 in [T3]) by the 
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change of variable y = {1 — z)/{l + z) ). The numerical 
evaluation of the integral gives approximately the value 
I2 — 11.2. This numerical evaluation results in the sur- 
vival probability whose asymptotical behaviour at t — ?> cxd 
follows: 



takes the inverse Laplace transform, 



exp(-x2^) 



and finally returns to L to get 



(36) 



where C = h/i^n^''^) « 0.64. 



Dt 



(30) 



Marginal distributions for the RWRW 

For the sake of completeness we give here the expres- 
sion for P(L, z), the joint probability density that a point 
on a track chosen at random falls in the interval of length 
L between its two crossing points with the x = const line 
and divides this interval into two parts whose length quo- 
tient is — z. This one is obtained by the variable 
transformation fromp(^i,Z2) = p{W\x)p{l2\x) with p{l\x) 
given by Eq. ([22]) and reads 



(1 



AL 



(31) 



The marginal distribution in z follows elementary |14) : 



p{z) 



r°° 1 

/ P{L,z)dL^- 
Jo 



1 



V^(l + z)' 



(32) 



This distribution corresponds to a very uneven distribu- 
tion of the starting point within the interval (the equidis- 
tribution would correspond to p{z) ~ (1 + ^)~^, and 
doesn't have a singularity at small z). This means, that 
in our case the particle starts close to one of the ends of 
the interval (the distributions of z and 1/ z are the same) , 
and therefore the FPT distribution decays faster than in 
the case of the equidistribution. 

To calculate the marginal distribution in i, one first 
passes to the distribution of the variable ^ = 1/L: 



1 + z 



exp 



(1 + Z)2 x"' 



takes the Laplace transform in ^ to obtain 



Piu,z) 



+ z) 



TT y/z{4:uz + a;2(l + z)2) ' 



(33) 



(34) 



performs integration in z (which is an elementary but 
tedious procedure) to obtain 



1 



x\l + z) 



-dz = 



^.{Auz + x'^il + z)^) \A7T^ 



(35) 



the Smirnov distribution, as anticipated. 

Note that neglecting the role of the initial position 
(which, as stated, is extremely unevenly distributed 
within the i-interval) would lead to a different result, 
namely to: 



27r3/2V^ 



1^ [^-K^Dt) 



(38) 



(see Eq. 2.3.15.5 in [IS])- Using the behavior of the 
modified Bessel functions for small values of its argument, 
^1/4(2) ~ ir(l/4)(z/2)-i/4 (see Eq. 9.6.9 in [IS]) we 
would get that $(t) cx r^l^ hke in the CTRW. 
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